Bulk black holes radiating in non-Z2 
brane-world spacetimes 

David Jennings^* Ian R. Vernon^j 
Anne-Christine Davis^l Carsten van de Bruck^§ 

^Department of Applied Mathematics and Theoretical Physics, 
CMS, University of Cambridge, Cambridge, CBS OWA, UK. 
Centre for Particle Theory, University of Durham, 
South Road, Durham, DHl 3LE, U.K. 
^Department of Applied Mathematics, 
University of Sheffield, Sheffield, S3 7RH, U.K. 

February 1, 2008 



Abstract 

In this paper we present a general asymmetric brane model involv- 
ing arbitrary energy transport to and from an embedded 4-D FRW 
universe. We derive a locally defined mass function for the 5D space- 
time and describe its time evolution on the brane. We then specialise 
our model to the two cases of graviton production in the early universe 
and radiating black holes in the bulk. 



*D. Jcnnings@damtp.cam.ac.uk 
^ I . R. Vernon@dur .ac.uk 
A. C.Davis@damtp.cam.ac.uk 
§ C . vandeBruck@sheffield .ac.uk 



1 



1 Introduction 



The idea of having extra, unseen, spatial dimensions to our universe has its 
origins with Kaluza and Klein [1,2] who proposed a five dimensional version 
of General Relativity that incorporated electromagnetism. Their idea was 
to include a small, compact, extra dimension. Our inability to observe this 
extra dimension was attributed to its scale being outside our experimental 
reach. The notion of small, hidden extra dimensions later reappeared with 
the advent of string theory and superstring theory where their existence was 
necessary for the finiteness of the theories. 

The modern extension of string theory is M-theory, where, in addition to 
one dimensional strings, we have higher dimensional objects called branes. 
Brane- worlds spring from a model suggested by Horava and Witten [3, 4] 
where the strong coupling limit of the £"8 x Eg heterotic string theory at 
low energy is described by eleven dimensional supergravity. An important 
new clement in this model is that one of the extra dimensions, the orb- 
ifold dimension, can actually be quite large. The reason for this is that the 
standard model particles are represented by open strings while gravitons are 
represented by closed strings, these open strings are constrained to finish 
on a brane while the closed strings have complete freedom. In brane-world 
scenarios [5-14] our universe is then viewed as a 3-brane embedded in a 
larger spacetime with one extra dimension which is not necessarily compact 
(see [15-17] for recent reviews). The standard model particles are attached to 
this brane, while the gravitons are free to move off it into the extra dimension, 
which is referred to as the bulk. 

The usual model consists of a bulk Einstein-Hilbert action together with 
a brane action. 



where V is the brane tension for our 3-brane universe and A is the bulk 
cosmological constant. As such the brane is viewed as a hypersurface in the 
5-D spacetime. Its induced metric and extrinsic curvature are defined as 



brane 



EH 




(1.1) 



Kab = h'^h^VcriD, 



(1.2) 
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where ua is a vector field normal to the hypersurface and capital indices run 
over the five dimensions. Since we assume that the standard model particles 
are confined to the brane, the matter contribution from our universe will take 
the form of a delta function in T^. Then, using Einstein's equations we may 
obtain the famous Darmois-Israel junction conditions [18] relating the jump 
in the extrinsic curvature across a hypersurface to its matter distribution on 
the hypersurface, 

[Kab] = —1^^ ^^A_B,branc " "^AB^ branc ^ • (1-3) 

The spatial components of the Israel-junction conditions for the brane pro- 
duce a modified FRW equation for the Hubble parameter [19-22], 

where we have assumed that the brane universe is homogeneous and isotropic 
and that a Z2-symmetry is present across the brane. This differs from the 
standard FRW equation in two ways: firstly we now have a term quadratic in 
the density, whose presence would dominate in the very early universe; sec- 
ondly there is a term which is referred to as the dark radiation and which 
would also become relevant in the early stages of the universe. The dark ra- 
diation term is actually the one remaining component of the projected Weyl 
tensor [23] for the bulk spacetime and so is an external feature determined by 
the content of the bulk, e.g. gravitons fiowing in the extra dimension. This 
extra term could provide a possible test for the above extra-dimensional sce- 
narios. Although many authors have assumed that C is constant, this will 
not be the case if there is energy fiowing to, form or across the brane. 

In this paper we develope a general formalism to describe the dynam- 
ics of an isotropic, homogeneous non-Z2 (or asymmetric) 3-brane in a bulk 
that possesses arbitrary energy transport to, from and across the brane. We 
then use these results to examine a variety of situations, including the case 
where the energy transport is due to the bulk black holes either side of 
the brane emitting Hawking radiation, in addition to the brane thermally 
emitting gravitons. The paper is organised as follows: in section |21 we sum- 
marise the important aspects of non-Z2 symmetric branes, before presenting 
the necessary formalism to determine their evolution in a general bulk. We 
derive the Friedmann equation along with the necessary equations for the 
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time-dependent dark radiation parameter C and non-Z2 parameter F, then 
we discuss the locally defined mass function, related in some cases to the 
black hole mass. In section El we apply our general equations to the cases 
of graviton emission by the brane and to Hawking radiation by both black 
holes, before presenting a numerical treatment and discussion of the various 
cases in section 0] Our conclusions are given in section El 



2 Non-Z2 branes 

In this section we develop the equations governing energy transport to and 
from a non-Z2 brane universe, where we do not identify the two bulks ei- 
ther side of the brane [24-27]. We then specialise to the case of radial null 
fiuids in the bulk. The presence of several overlapping null fluids in the 
bulk makes the overall description problematic and an analysis using Vaidya 
metrics [28-30] will not work. The natural coordinate system to use would 
seem to be a double null system of ingoing and outgoing null rays, but this 
presents several technical difficulties which makes analysis awkward. It is 
convenient for our purposes to began with the Einstein equations in 5-D, 
using a coordinate system based on the brane - as opposed to tackling issues 
of extrinsic curvature in a bulk-based system. 

2.1 General description 

We work with a five dimensional metric describing a standard 3-brane with 
the usual homogeniety and isotropy for its spatial dimensions, 

^ -n'^{t, z)dt^ + a^(t, z)-fijdx'dx^ + z)dz^, (2.1) 

and a total 5D energy-momentum tensor which we decompose as, 

branc, vac vac branc, matter matter ' \ ' / 

We shall assume separate cosmological constants on the brane and in the 
bulk, and a perfect fiuid on the brane. At present we shall make no assump- 
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tions concerning the energy-momentum tensor in the bulk. 

7^i,..ane,vac = ^-^diag{-V, -V, -V, -V, 0) , 

= H-p,o-r-^ 

B,bulk, vac ± CUagl 2' ^2' ^2' ^2' ^2''' 

A/ My A/ My A/ 

^B.^^ane, .aue. = ^diag( -p, p, p, p, 0) . (2.3) 

Here the ± refer to the two bulk spacetimes either side of the brane and 
we shall work without the usual Z2-symmetry across the brane. To aid the 
discussion we define, for any quantity X, its jump and average across the 
brane, 

AX = Xq+ — Xq- , 
X=(Xo++Xo-)/2, (2.4) 

The Einstein equations 

Gab = k^Tab, (2.5) 

have Gab given by [19,20], 

^ ^ I 0- I 0- b\ f a" a! ( o! 

Goo = 3 - - + - -- - + - 
\a\a I o^\a a 

^ f a' f a' ^n'\ h' f n' a'\ ,„ , 

Gij = 1^7., - - + 2- -- - + 2- +2- + — + 

ah „ a ^ab hb b\ 
__2--2- + -- - - kj,„ 
an a ab nb b 



05 





- 3(^(^ + ^l-f,(^+«f«-?)l-4). (2.6) 



a a \ a nil cP- 



2.2 Active source equation - The jump conditions 

We may use the fact that only the second derivatives with respect to z contain 
delta functions and integrate the (OO)-equation and the (ij)-equation over the 
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small interval [— e, +e] , which contains the brane at z = 0. This provides us 
with the jump conditions across the brane in terms of the matter content on 
the brane [19], 

An' = —nobo(3p+2p-V), 

Aa' = -'^aMP + V), (2.7) 

where the primes denote derivatives with respect to the transverse coordinate 
z, and Xq is the quantity X evaluated on the brane ai z = 0. For convenience 
we go to the temporal gaussian coordinates where, 

bo = no = 1, 

no = 0. (2.8) 

For the sake of clarity, we now drop the subscript 0. In addition, all unspec- 
ified terms Tab shall henceforth refer to the bulk content. 



2.3 Passive source equation 

The Darmois- Israel junction conditions on their own are incomplete when the 
Z2-symmetry across the brane is removed. It is thus necessary to supplement 
them with a passive gravitational source equation describing the motion of 
the brane within the background spacetime. The origin of this passive source 
is found in the divergence equation for the brane universe [31], 

VaTII = (2.9) 

where is a force density localised on the brane. We then consider the 
orthogonal component of this field, 

/ = riBf^. (2.10) 

After some tensor manipulation, the dynamical equation is integrated across 
the brane and we obtain the passive gravitational source equation, 

KLKab = 7- (2.11) 

It has been observed [31] that this is analogous to Newton's second law of 
motion, with Tab taking the place of mass and Kab being the acceleration. 
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In our current setting we may obtain the explicit form of the passive 
equation by analysing the (55)-equation. We find that it takes the form, 

^(p +V)- 3-{p -V) = ^- AT55, (2.12) 
which for a Z2-symmetric brane is satisfied trivially. 

2.4 Transverse equations 

The (05)-equation provides us with the equation governing energy flow off or 
onto the brane. With the help of the jump conditions, which relate transverse 
derivatives to the energy content of the brane, we find that this takes the 
form, 

p + 3H{p + p) = AT05, 

H ^ -. (2.13) 

This form is quite natural, with the usual energy conservation replaced with 
one containing a 'leaking' term, AT05, which corresponds to the net energy 
flux into the extra dimension. 

The average of the (05)-equation tells us how the geometrical asymmetry 
across the brane evolves with energy transfer from one side to the other. We 
find that, 

-dt^ = Ti/H-'^To,, (2.14) 
a 6 

which will be needed later. 

2.5 The FRW equation for the brane 

Prom the (55)-equation we may now arrive at an acceleration equation for 
the brane scale factor, 

H' + l + ^ = ^i-(^) [^(3p-p) + p(3p + p)] + e-^T55, 




(2.15) 
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Where B is the geometrical asymmetry contribution. We now define the 
auxihary field which represents a generalised dark radiation term, via the 
FRW equation, 

= /3p2 + 27p-4 + X + A, (2.16) 



a 



where /3,7, A are constants. Using this together with equation ()2.15j) we find 
that the auxiliary field dynamically obeys, 

X + ^Hx + 2h(^T,, - 0) = -^(p + V)Aro5. (2.17) 
o io 

From its form we see, in the absence of asymmetry, that it behaves like 
radiation coupled to matter on the brane. Using the acceleration equation 
for the brane scale factor and the defining equation for x we require that the 
constants in the definition of x have the following values, 

/3 = kV36, 
7 = 

b 

It is now possible to split this generalised dark radiation term, x, into two 
natural parts. To describe the asymmetric geometry across the brane we 
define the function F{t) as, 

a'= ^ 3 ■ 2.19 

(p + Vja'^ 

We now decompose x iiito two parts and in the process introduce the Weyl 
term C{t) which is the usual dark radiation factor in the FRW equation. 

This gives the FRW equation for a non Z2-symmetric brane in a completely 
general bulk. 



(p + Vya^ 



= Pp' + 2^p-^ + ^+ ; fiXr; ' + A, (2.21) 
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2.6 Evolution of parameters 

F{t) represents the asymmetry of energy across our brane boundary, and 
consequently we expect its development to be decided by transverse flows to 
the brane. With this in mind we return to equation ()2.14|) and find that, 



+ -^AAa') - Ha' AT,, - ^a\p + V)T,,. (2.22) 



It is then a simple matter to deduce that the Weyl term C{t) must obey, 

= + ^AAa^) - - ^a\p + V)AT,,. (2.23) 

These equations, together with, 

H - Pp +2^p-- + -^+ ^p^yy^, +A, 
p = -3iJ(p + p) + ATo5, (2.24) 

constitute a complete set of equations describing the cosmological evolution 

of the brane-world model. 

We now make the observation that, 



55,- 



and so the two quantities -^C + |F and -^C — |F are intrinsic in some 
way to the + and — bulks respectively. We shall demonstrate that these two 
quantities correspond to the energy content of their respective bulks, first by 
considering Vaidya spacetimes and then in a more general setting. 

2.7 Vaidya spacetimes 

The Vaidya spacetimes [32, 33] generahse AdS-Schwarzschild and describe 
the inflow or outflow of radial null fluids. The metric for such a spacetime 
takes the form, 



ds^ = —f{v, r)dv'^ + 2edvdr + r'^jijdx^dx^ 
A 2 2m{v 
6 r 



/ = -ii.^-^, (2.25) 
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where v — constant describe radial null rays which are ingoing for e — 1 and 
outgoing for e = — 1. Wc may derive the FRW equation for a brane in this 
spacetime by using the 5-D junction conditions [18], 



[Kab] = i<^{Tab - ^THab), 



(2.26) 



where Hab is the induced metric on the brane. If we do not impose Z2- 
symmetry we obtain a FRW equation [30] , 



k 2m (^Am + ^AAa^)^ 



(p + y)2a8 



+ A, 



(2.27) 



and so we may infer that, 



C 



2m, 



K 



—F = Am. 
3 



(2.28) 



We then write the equations governing the evolution of the cosmological 
parameters as. 



k 2m{t) (12Am + AAay 



36''''' ' ' ' a? of 



+ 



16K4(p + T/)2a8 ' 



Am 



m 



Am + 



p + V 



To5,± 3""^" ^55,±, 



^(Am + ^a^) - ^Ha^AT,, - |a^(p + V)To,, 
p+V 12 6 y 

^(Am + ^a*) - ^Ha'T,, - ^a\p + V)AT,,, 
p+V 12 6 3d 



where we have imposed the Randall- Sundrum tuning conditions [13, 14], 

0, (2.29) 



A + 



6 



which means the cosmological constant on the brane vanishes. 
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2.8 Locally defined mass function 

We now demonstrate that this identification of the bulk mass parameter with 
linear combinations of C and F may be extended to the case of a general 
bulk with the usual homogeneous and isotropic brane boundary. If we begin 
with the metric (j2.H) and introduce the following Kodama vector field [34], 

A;^ = (^,0,0,0,4), (2.30) 
no no 

we find that this vector field is conserved: 

v^^i = (v^fc^),A = 0. (2.31) 

We now use the fact that T^.^ = to define the following conserved current 
in the bulk spacetime, 



= ^^T^K''. (2.32) 

It is easily checked that OaJ^ = and so we may define a conserved charge 
M{t, z) = f J^dz satisfying, 

M(t,z) = -J^{t,z), (2.33) 

where J^{t, z) is the flux of M out through the z =constant hypersurface at 
the point (t, z). Explicitly, in terms of our metric components, this becomes, 

M{t,z) = -{Tf-^-Tf^)a{t,zY. (2.34) 
a a 

This expression is now evaluated on the brane from either side. For example, 

M(t,z = 0+) = -(Tl^H{t)-T^^^)a\ (2.35) 

We then make use of the jump conditons and the definition of the asymmetry 
function F, 

K-«-)/2 = -^{P+V), 
11 



to arrive at the result, 

M{t, . = 0^) = [^^-^^^ - Y^\P + V)^ To5,+ - Ha'n,,^. 
Similarly, approaching the brane from the other side we obtain, 

M{t,z = 0_) = [^^-^^P^ + Y^'iP+V?j To5,_ - Ha'T,,,_. 

We then conclude that the mass function defined locally at each point in the 
5D spacetime obeys, 

M(i,z = 0±) = ^C±h, (2.37) 

which provides us with the physical relation the two components of x have 
with the bulk spacetimes. 



3 Specific bulks 

We may now specialise our dynamical equations to different bulk contents. In 

order to describe the system completely we need to specify four independent 
physical quantities, Tq^^±, T55 j-, to describe the effects the two different bulks 
have on the brane. These are then combined together into AT05, AT55, T05, 

T55- 

If we consider a null fluid with energy momentum tensor given by. 

Tab = crkAks, (3.1) 

where kA is a radial null 5- vector in the bulk then wc may take to be the 
tangent vector to the brane and normalise so that TabU^u^ = o") where a is 
the out- flux of gravitons as measured at the brane. With this we find that, 

= a. (3.2) 
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3.1 Brane radiating gravitons 



This model has been considered recently by several authors in both the sym- 
metric [28, 29, 35, 36] and asymmetric [30, 37] cases. For this spacetime we 
have that (T± = ±^K^p^, where alpha is a dimensionless constant, and the 
evolution of the mass parameters is governed by. 



m± 



11 



K 



p + V 



-Ha'' 



a 
12 



2 2 
K p , 



(3.3) 



where we have used. 



T55 



^22 

6 ^ ' 
AT. 



55 

" 2 2 
12'^^- 



0, 



The equations for the full model become, 

2m(t) 



//2 

P 

Am 
m 



-3i/(p + p)- 



+ 



-(AmH a^), 

p + K ^ 12 ^' 



4 2 4 2 

an p a K 



36 



6 



{p^V)-H). 



4\2 



;i2Am + AAa ^ 



(3.4) 



(3.5) 



These equations have been studied in [30], and in [37] where the global nature 
of the bulk spacetimes is also examined. 



3.2 Hawking radiation from bulk black holes 

We now consider the situation of an asymmetric brane model with a black 
hole in the two bulk spacetimes either side of the brane. We take the black 
hole mass parameters to be different from one another and include a possible 
difference in the bulk cosmological constants. The difference in mass param- 
eters provides a contribution to the FRW equation, whose presence would be 
strong in the early universe. It has previously been assumed that these mass 
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parameters remain constant throughout the model's history, but the black 
holes can radiate thermally through Hawking radiation and vary the dark 
radiation and asymmetry terms in time. Since we assume standard model 
particles are constrained to the brane, and the bulk is free of any scalar fields, 
the only channel remaining for Hawking radiation is through bulk gravitons. 

We now calculate the evaporation rate for these bulk black holes. We as- 
sume that a non-zero cosmological constant in the bulk does not substantially 
affect the evaporation rate. Following [38] we make a spherically symmetric 
approximation for the black holes and consider the metric, 

ds^ = -f{r)dt^ + J^ + r^dn^ (3.6) 

with, 

m = 1-^' (3-7) 

For us To = 2m, which may be related to the ADM mass, M via (see [39]), 

ro = 2m = (3.8) 

In D dimensions, the rate of decrease of a radiating black hole in the high 
frequency limit is given by [40], 

-j^ = -goaDAoT'^, (3.9) 

where go is the number of degrees of freedom, a^, is a generalised Stefan's 
constant. Ad is the geometrical area of the event horizon and T is the tem- 
perature of the black hole as measured at infinity. In 5 dimensions with only 
graviton emission, this becomes 

f^-^C.K'.^. (3.0) 

We derive the black hole temperature by considering the periodicity of the 
Euclidean time for the metric. The Hawking temperature of the black hole 
is then given by, 

= (3.11) 
47r 2nro 
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Consequently we find that the rate of decrease of the mass parameter near 
infinity is given by, 

fdm\ ( 80C(5)k^ \ 

To arrive at the flux of gravitons at the brane we consider, as in section 2.7, 
a radially outgoing null fluid with the Vaidya metric, 

ds'^ = -fdv^ - 2dvdr + r^dn^. (3.13) 

The energy momentum in the bulk is again of the form Tab — crkAks which 
is normalised as, 

kAu"^ = 1, (3.14) 

so that a becomes the flux of gravitons measured on the brane. From the 
Einstein equations we flnd that, 

dm K\-a ^3_^^^ 



dr ?>v 



where v — and 



fy^^f^ + f- f. (3.16) 



If we now make an identification of these two systems near spatial infinity 
then we assume that the rate of decrease of the mass function at the brane in 
the Vaidya spacetime is equal to the derived rate of decrease of an evaporating 
black hole in AdS-Schwarzschild. Consequently we obtain, after replacing r 
with a for the scale factor. 



^2 f dm ^ 



-a — 



3 a^f 



2A0Ci5){y/dFTf - d) 

a3/(67r2)2(27r)5m ' ^ ^ 

which provides us with the energy fiux at the brane from the radiating black 
hole. We note that this evaporation rate contains a dependence on the brane's 
motion through the spacetime corresponding to a red-shift of the graviton 
fiow as measured on the brane. 
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We may now deal with the model consisting of two black holes, one in each 
bulk region, and combine together their respective fluxes { a± ) ai the brane. 
For this situation we can make the further assumption that the outflows and 
inflows do not interact too strongly with each other and we may analyse the 
ingoing fluxes and the outgoing fluxes separately. We may then write the 
total energy flux either side of the brane as, 

T5+ = a_ - a+, 

Tl = a.-a+. (3.18) 

Consequently we have that, 

AT" = 0, 

% = Tl^_ = Tl = a_-a+, (3.19) 

which merely states that the gravitons do not accumulate on the brane. 
The pressure terms for this system are given by, 

ATss = 0, 

T55 = |(7_ -(T+|. (3.20) 
We would then obtain the following set of equations, 

rr2 '^^ . oT.N k 2m(t) (12Am + AAa^)2 

p = -3Hip + p), 

.2\ 2 



Am = (^^) a'{p + V)T^, 



ro A A k2 

m = —^{Am + —a'')-—Ha^n5. (3.21) 

We note that the process of radiating black holes tends to increase the mass 
asymmetry across the brane, while the average mass may increase or decrease, 
depending on the relative sizes of the brane's expansion rate and the bulk 
asymmetry. 

A novel possibility is a model where there is initially a pseudo Z2-symmetry 
present across the brane in that both bulks have the same cosmological con- 
stant and mass parameter. However, if only one side contained an actual 
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black hole, then the resultant radiation would naturally create an asymme- 
try by its flux of gravitons. 

We now pass to dimensionless variables, where we note that the variables 
possess different dimensions depending on whether is or ±1. Using X to 
signify dimensionless equivalent of X, our equations become. 



dp 

di 
(iAm 

di 

dm 
di 



k 2m (Am + 3AAa^)2 



4(p + lyd^ 



-AHp, 

24f^d\p+l), 

6(Am + 36AAa^) 



P+1 



l2Ha^t 



55- 



(3.22) 



It is also straightforward to include the terms governing graviton production 
on the brane. The resultant system possesses certain simplifications such 
as the fiow of gravitons being purely radial and the assumption that once 
energy enters a given bulk it immediately increases the black hole mass. We 
provide simulation results for this model in the next section. 



4 Numerical analysis 

In this section we present the results from numerically simulating bulk black 
holes that radiate. We considered the four models obtained by turning off 
or on the brane flux or the black hole fluxes. It is then possible to deter- 
mine the separate contributions of each process to the overall cosmological 
model of a radiating brane moving through a bulk containing radiating black 
holes. In what follows, we discuss the black hole features separately from the 
cosmological variables on the brane universe. 

4.1 Black hole flux and mass 

We solved numerically equations ()3.22|) with a range of initial black hole 
masses, m+ and m_. Each simulation was then repeated with a difference in 
cosmological constants between the two bulk spacetimes. 
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In figure 1 we show the black hole fluxes as measured at the brane for 
an average black hole mass of m = 6 and a difference of Am = 2. The 
solid curves (cr-t) are for the model containing both Hawking radiation in the 
bulk and graviton production on the brane. The dotted curves {(t^^±) are 
the results without graviton production. We note that the fluxes are only 
on when the brane is outside the event horizon. The resultant motion of the 
brane causes a reduction in the energy flux of gravitons impinging on the 
brane, as it moves further and further away from the black holes. The fluxes 
(7+ and (T_ gradually become similar in magnitude and so the variation in 
the dark radiation and asymmetry terms is substantially reduced. 



M=0,k=0,Am=2,m=6,V=0/600,a=0/2.5 




t 



Figure 1: Black hole fluxes according to a brane based observer, initially 
with Am = 2 and m = 6. The solid lines {(t±) represent the fluxes on the 
right/left side of the brane for the black hole radiating plus brane emitting 
gravitons case. The dotted lines represent the black holes' fluxes when 

the brane does not emit gravitons. 

For a larger mass difference across the brane, it is easier for the smaller 
black hole to evaporate completely and the fluxes of the two black holes do 
not cancel, as is shown in figure 2. The flux of the smaller black hole initially 
decreases to a minimum and then diverges as it evaporates. We also note 
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that the addition of graviton production on the brane enhances black hole 
flux in the red-shift phase and diminishes it during the evaporation burst. 
This effect occurs due to the brane's path, a(t), being perturbed by the 
gravitons leaking out of the brane. The addition of a small difference in 
the cosmological constants across the brane results in a slowing down of the 
previously mentioned behaviours. 




t 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Figure 2: Black hole fluxes according to a brane based observer, initially with 
Am = 4.5 and m = 6. Notation is the same as figure 1. 

The evolution of the average black hole mass, in the brane's frame, is 
not monotonic in the scenario of pure black hole radiation, as there are 
competing asymmetry and Hubble terms. Initially, the expansion of the 
universe dominates over the mass asymmetry and the average black hole 
mass decreases. From equation ()2.29|) the average mass becomes stationary 
if, 

3Am+^a^ = (/t^o + ^a^)^, (4.1) 

after which the average mass will increase. The addition of graviton produc- 
tion in the brane does not alter this trend to any great extent, see figure 3. 
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The addition of a small, positive AA across the brane tends to strongly in- 
crease the average black hole mass in the radiating black hole model, but has 
a minor effect in graviton production model, figure 4. 




Figure 3: The behaviour of m{t) in the four models. The black solid line 
represents the case where both the brane and the black holes are radiating, 
the green/red line the case where only the brane/black holes are radiating, 
and the dashed line where no radiation occurs 

In the case of pure black hole emission the mass asymmetry always in- 
creases in time, as can be seen from equation ()2.29j) this competes against the 
emission of gravitons from the brane which tends to reduce the asymmetry. 
Because the flux of gravitons out of the brane depends quadratically on the 
density of the brane, it can initially dominate over the black hole radiation, 
but later has little or no affect on the development of the mass asymmetry, 
figure 5. If we insert a positive difference of bulk cosmological constants in 
the black hole radiation model the increase in asymmetry is unchanged, as 
opposed to the radiating brane model where the asymmetry reduces more 
rapidly, figure 6. 
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Figure 4: The behaviour of m{t) in the four models with AA = 0.1. Labelhng 
is the same as figure 3. 

4.2 Cosmology on the brane 

We now turn to the effect these radiating black holes have on the cosmo- 
logical parameters in the brane. Their presence manifests itself primarily 
through the dark radiation and asymmetry terms, and indirectly through 
the evolution of p{t) and a{t). 

In figure 7 we have plotted the fractional contribution of the dark radi- 
ation term to the FRW equation. We consider its variation over the four 
models obtained by switching on or off the radiating black holes or the radi- 
ating brane. We see that having radiating black holes in the bulk tends to 
suppress the dark radiation as opposed to the radiating brane, that increases 
the dark radiation by transfering energy to it from graviton production. In 
contrast, radiating black holes tend to magnify the contribution of the asym- 
metry term as opposed to brane radiation which has little affect, figure 8. 

The addition of a difference in cosmological constants tends to provide an 
effective late time cosmological constant in the FRW in the asymmetry term, 
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Figure 5: The behaviour of Am(t) in the four models. Labelhng is the same 
as figure 3. 

and is thus constrained experimentally to be small. It also has the effect of 
increasing the contribution of dark radiation over time, by enhancing the two 
graviton processes present. 

The scale factor is modified in different ways by the brane and black hole 
radiation. In general the black hole radiation tends to give a larger scale 
factor for the universe, while radiation from the brane tends to reduce it. 
In the presence of a small cosmological constant we may compare the four 
models in terms of their scale factors, as in figure 9. The Hubble parameter 
is decreased by both radiation processes, see figure 10, with brane radia- 
tion dominating the suppression at early times until the black hole radiation 
increases in strength. 



5 Conclusion 

In this paper we developed the general framework necessary to describe an 
asymmetric brane universe with arbitrary bulk content. We highlighted the 
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M=0.1,k=0,Alii/m=0.5,V=0/600,a=0/1.5 
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Figure 6: The behaviour of Am{t) in the four models with AA = 0.1. La- 
belhng is the same as figure 3. 

fact that in the absence of Z2-symmetry, it is necessary to include a passive 
source equation, describing the motion of the brane within the background 
spacetime. Through the FRW equation we were able to introduce a gener- 
alised dark radiation term, which we showed was composed of two natural 
pieces. One piece is the familiar dark radiation term, while the other con- 
tained the asymmetry terms for the model. Their evolution in time indicate 
that the expressions ± are intrinsic to the two different bulks and a 
calculation in a Vaidya spacetime revealed that they correspond to the mass 
parameters. We then arrived at a locally defined mass parameter which gen- 
eralised this identification with the dark radiation and asymmetry terms in 
the FRW equation. 

In section 3 we retricted our equations to two cases of particular interest: 
graviton production on the brane and radiating bulk black holes. We calcu- 
lated the flux of gravitons from these bulk black holes as measured on the 
brane and performed numerical simulations of the system. It was found that 
radiating black holes tend to increase the asymmetry term in time, while the 
brane radiation reduces it. The effect on the dark radiation was less straight- 
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Figure 7: The fractional contribution of dark radiation to H^. Labelling is 
the same as figure 3. 

forward with a competition between cosmological expansion and asymmetry 
terms. 

The two radiation processes produce corrections to the scale factor and 
the Hubble parameter (figures 9 and 10), while the black hole radiation 
amplifies the asymmetry term in the FRW equation. The analysis we have 
presented can be used to determine the values of m and Am at the time 
of nucleosynthesis, and thereby to impose strict constraints on the initial 
conditions of the system, but we leave this to future work. Other interesting 
investigations in this area could include an analysis of non-radial graviton 
flows and the effects associated with having scalar fields present in the bulk. 
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